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We calculate the neutrino nucleus cross sections for charged lepton production relevant for the 
experiments proposed with the stopped muon neutrinos using neutron spallation source facility. 
The calculations are done in local density approximation taking into account Pauli blocking, Fermi 
motion effects and renormalization of weak transition strengths in the nuclear medium. The effect 
of Coulomb distortion of the lepton produced in charge current reactions is taken into account by 
using the Fermi function as well as in a model where an effective momentum has been used for the 
lepton moving in the local Coulomb field of the final nucleus. The numerical results for the neutrino 
nucleus total cross sections averaged over Michel spectrum are presented for various nuclei. 
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I. INTRODUCTION 

The study of neutrino nucleus cross sections at low energies is important in understanding various physical and 
astrophysical processes. The low energy neutrino beams (E v < 52.8MeV) are generally obtained from the muons 
decaying at rest. The neutrino energy spectrum obtained from the muons decaying at rest is given by 



<t>{E v . ) = ^zE 2 Ve (E -E Vb ), E = 52.8MeV 



(1) 



This is known as Michel spectrum and is shown in Fig.l. 




FIG. 1: Michel spectrum for v e obtained from muon decay at rest. 



This neutrino energy spectrum and its energy range is similar to the energy spectrum and energy range of neutrinos 
coming from the core collapse supernova 1]. This similarity in the energy range and spectrum of the supernova 
neutrinos with the muons decaying at rest opens up the possibility of connecting the ground based neutrino nuclear 
experiments with the study of neutrino nuclear cross sections in supernova. Such a study will also be useful in 
understanding the r-process nucleosynthesis leading to the formation of heavy elements in the interstellar medium^]. 

The low energy neutrino beams can also be obtained at the Spallation Neutron Source(SNS) facilities where low 
energy neutrino nucleus reaction experiments can be performed 3] . At these facilities the accelerated protons hit a 
nuclear target and produce tt~ and ir + particles which give rise to muon and neutrinos as their decay product. There 
is a small flux of D e as a result of [i~ decay as most of pT are captured by nuclei. On the other hand fi + are not 
captured but are stopped in nuclei and decay at rest giving rise to substantial flux of v e . These neutrinos can be 
used to study neutrino nuclear cross sections. Experiments to study the feasibility of making such neutrino- nucleus 
cross section measurements using many nuclear targets are proposed at Oak Ridge National Laboratory(ORNL) using 
stopped muon neutrino facility at the Spallation Neutron Source(SNS). Such a facility is expected to produce ss 10 15 
neutrinos per second of each flavor of v e , za, and with mono-energetic spectra for Vn and continuous spectra for 
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v e and Vp. Out of these flavors of neutrinos, v e produces electrons which are observed to study the neutrino nuclear 
cross sections for charged current reactions. 

The electron production from nuclei induced by low energy v e beams has been studied in past by many authors. 
A review of theoretical calculations and experimental results have been given by Donnelly and PecceiQ, Kubodera 
et al.0 and more recently by Kolbe et al.p. It is proposed that these calculations are done for a large range of 
medium and heavy nuclei and attempts be made to study them experimentally. Such a list of nuclei where the study 
of neutrino nuclear cross section is of some interest has been recently prepared by Avignone and Efrcmcnko 7] . 

In this paper we have studied neutrino nuclear cross sections induced by low energy v e obtained from muons 
decaying at rest for the nuclei proposed by Avignone and Efremenko0. The calculations have been done in local 
density approximation taking into account Pauli blocking, Fermi motion effects and renormalization of weak transition 
strengths in the nuclear medium. The method has been successfully applied to study the various electromagnetic 
and weak processes in nuclei at low and intermediate energies|8j-|l4|. The effect of Coulomb distortion of the lepton 
produced in charge current reactions is taken into account by using the Fermi function F(Z,E e ), where Z is the 
atomic number and E e is the outgoing electron energy as well as in a Modified Momentum Approximation(MEMA), 
where the effect of Coulomb distortion is incorporated by modifying the momentum and energy of charged lepton in 
the Coulomb potential of the final nucleus 0,0. 

In section-II, we describe the formalism. In section-Ill, we present the numerical results for the total cross section 
(j(E) vs the neutrino energy E and the flux averaged cross section and discuss the nuclear medium effects in some 
detail. A brief summary and conclusion of our results is given in section-IV. 



II. FORMALISM 

In this section we derive an expression for the total scattering cross section <j(E v ) for the charged current(CC) 
reaction 



X 



a Y 



where Zi{Zf) is the charge of initial(final) nucleus. 

The basic I'e-neutron reaction taking place in %.X nucleus is 



(2) 



(3) 



where k and p are the four momenta of the incoming neutrino and neutron and k' and p' are the four momenta of 
the outgoing electron and proton respectively. 

The matrix element for the basic neutrino process on free nucleon (Eqn.3) is written as 



where 



v2 



= u(p')[F? '(gV 



(4) 



(5) 



The form factors F^(q 2 ), F^ (q 2 ) and F\(q 2 ), where q 2 is the four momentum transfer square i.e. (fc — fc') 2 , are 
isovector form factors and are written as 

F?{q 2 ) = Fl(q 2 ) - F?(q% F? \q 2 ) = F^q 2 ) - F?(q% 
FX(q 2 ) = F A (q 2 ) 

where 



1 - 



q 



AM 2 
q 2 



AM 2 



G p E n (q 2 )-^G^(q 2 ) 

[G p M n {q 2 )-GT{q 2 )] 



(6) 
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(i-ife) 2 



G P M (q 2 ) = (1 + P P )G p E (q 2 ), G n M {q 2 ) = p n G p E (q 2 ); 



„2 



GW) = (^)»nG E (q 2 )Zn; & 



s 4M 2 ' r " ^ s " i_a -a 



li\/ 2 



= 1.792847, \i n = -1.913043, M v = OMGeV, and A„ = 5.6. 
The isovector axial vector form factor Fa{Q 2 ) is given by 

where M A = IMGeV; F A (0)=-1.26 

The double differential cross section <r a (E e , \k'\) for the basic reaction described in Eqn.(3) is then written as 

^ 2 

a (E e , \k'\) = G F 2 cos 2 9 ]k ' 1 *i n *f p XX\T\ 2 6[q + E n E p ] (7) 

where ££|T| 2 is the square modulus of the transition amplitude given by 

££|T| 2 = L^J^ with 

'g^u + itnu\ a k x k' a ] and 

= J2 j ^ r ( 8 ) 

(70(90 = E Uc — E e ) is the energy transferred to the nucleon. 

In a nucleus, the neutrino scatters from a neutron moving in the finite nucleus of neutron density Pn{r), with a 
local occupation number n„(p, r). In the local density approximation the scattering cross section is written as 

a(E e ,\k'\) = J 2drdpj^n n (p,r)a (E e ,k') (9) 

where ao(E e , \ k'\) is given by Eqn.(7). The neutron energy E n and proton energy E p are replaced by -E^Op]) and 
E p (\p + q\) where p is now the momentum of the target neutron inside the nucleus. In the nucleus the neutrons 
and protons are not free and their momenta are constrained to satisfy the Pauli principle, i.e., p n < pf„ and p' p (= 
\Pn + q|) > Pf p , where pp n and pp p are the local Fermi momenta of neutrons and protons at the interaction point in 

the nucleus and are given by pp n — [37r 2 p„(r)] 3 and pp p — [37r 2 p p (r)] 3 , p„(r) and p p (r) are the neutron and proton 
nuclear densities. Furthermore, in nuclei the threshold value of the reaction i.e. the Q-value of the reaction has to be 
taken into account. 

To incorporate these modifications, the S function in Eqn.(7) i.e. S[q + E n — E p ] is modified to 8[q a + E n (p) — 
E p {p + q) — Q] and the factor 



/ 



j0r Mp,r) ^; s[qo+En ~ Ep] (10) 



occurring in Eqn.(9) is replaced by — (l/Tr)lmUN(qo, q), where UN(qo,q) is the Lindhard function corresponding to 
the particle hole(ph) excitation shown in Fig. (2) and is given by 



dp M n M p n n (p) [1 - n p (p + q)} 
u N {qo,q) I {27V )3 En E p q + E n (p) - E p (p + q) + ie [ > 



where qo=E Ve - E e - Q. 
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FIG. 2: Diagrammatic representation of the neutrino self-energy diagram corresponding to the ph-excitation leading to v e + n — > 
+p in nuclei. In the large mass limit of the lV~B(i.e. Mw — * oo) the diagram 1(a) is reduced to 1(b) which is used to calculate 

1 1 2 



\Ty in Eqn.(8). 



The imaginary part of the Lindhard function is then derived to be: 

ImU N (q ,q) 



1 MpMn [p ... 



(12) 



q 2 < ,i?F 2 — Qo < E Fl and 



< E Fl , where = \/W„ 2 + M n , £ 



^2 



p Fp 2 + M p 2 and 



Max 



M n , E F2 - q Q , 



-9o + 



The threshold value, Q, for the neutrino reaction is generally taken to be the Q value corresponding to the lowest 
allowed Fermi or Gamow Teller transitions. However, in some cases the Q value corresponding to the ground state to 
ground state(gs-gs) transition is also takenflfij. 

With inclusion of these nuclear effects the cross section a(E u ) is written as 



o{E v ) = - 



2G F 2 cos 2 C 



r 2 dr 



p e 2 dp e / d(cos9) 



1 



E V E, 



■L^J^ImU N [E Ve -E e -Q, q]. 



(13) 



In the nucleus the strength of the electroweak coupling may change from their free nucleon values due to the presence 
of strongly interacting nucleons. Conservation of Vector Current (CVC) forbids any change in the charge coupling 
while magnetic and axial vector couplings are likely to change from their free nucleon values. These changes are 
calculated by considering the interaction of ph excitations in the nuclear medium in Random Phase Approximation 
(RPA) as shown in Fig. 3. The diagram shown in Fig. 3 simulates the effects of the strongly interacting nuclear medium 
at the weak vertex. The ph-ph interaction is shown by the wavy line in Fig. 3 and is described by the n and p exchanges 
modulated by the effect of short range correlations. 

The weak nucleon current described by Eq.(5) gives, in nonrelativistic limit, terms like FaBt^ and iF2^jj-r + which 
generate spin-isospin transitions in nuclei. While the term iF2^jT + couples to the transverse excitations, the term 
Fa(jt + couples to the transverse as well as longitudinal channels. These channels produce different RPA responses 
in the longitudinal and transverse channels when the diagrams of Fig. 3 are summed over. This is illustrated by 
considering a term like Fact' in Eq.(5). One of the contributions of this term to the hadronic tensor J 4 - 7 in the 
medium is proportional to F^Sijlml/N which is now split between the longitudinal and transverse components as 

FlSijlmUN -> F\ [q ; qj + (% - q s qj)] ImU N (14) 
The RPA response of this term after summing the higher order diagrams like Fig. 3 is modified and is given by Jr P a'- 



J ij -> J l l PA = F 2 A ImU N 



qiqj 



1 - U N Vi 1 - U N V t 



(15) 
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FIG. 3: Many body Feynman diagrams (drawn in the limit Mw — * oo) accounting for the medium polarization effects 
contributing to the process z/ e + n — > e~ + p transitions 



where Vi and Vt are the longitudinal and transverse part of the nucleon-nucleon potential calculated with n and p 
exchanges and are given by 



f 2 
mi 



Vt(q) 



'A 2 



A 2 -g 2 



Ap 2 ~ 
A 2 - « 2 



(16) 



Ajr = 1.3Gey, C p = 2, A p = 2.5GeT^, and m p are the pion and p masses, and g' is the Landau- Migdal parameter 
taken to be 0.7 which has been used quite successfully to explain many electromagnetic and weak processes in 

nuclei^, 

This modified tensor Jr P a when contracted with the leptonic tensor Lij gives the contribution of the F\ term to 
the RPA response. The contribution of the time component of the hadronic tensors like J 0i and J 00 arc of higher- 
order in and are not important at low energies considered in this paper. 

The effect of the A degrees of freedom in the nuclear medium is included in the calculation of the RPA response 
by considering the effect of ph-Ah and Ah- Ah excitations as shown in Fig. 3(b). This is done by replacing Un by 
Un — Un + Ua, where Ua is the Lindhard function for Ah excitation in the medium and the expressions for Un 
and Ua are taken from The different couplings of N and A are incorporated in Un and Ua and then the same 
interaction strengths V/ and V t are used to calculate the RPA response. This is discussed in some detail in Ref and 
more recently in Ref.[ll| by Nieves et al. 

One of the important aspects of charge current neutrino interactions is the treatment of Coulomb distortion of the 
produced lepton in the Coulomb field of the final nucleus. At low energies of the electron relevant to j3 decays in 
nuclei the Coulomb distortion of electron in the nuclear field is taken into account by multiplying the momentum 
distribution of the electron in Eqn.(13) by a Fermi function F(Z, E e ), where F(Z, E e ) is given by[ 



F(Z,E e ) = 



1-3(1-70) 



f(Z,E e ), 



where 



f(Z,E e ) =2(l + 70 ))(2p eJ R) 



-2(1 



- 70 ) |r( 7 Q + »7)| 2 

(r(2 7o + 1)) 2 



Here R is the nuclear radius and 70 = yl — (aZ) 2 , r\ — stae. This approximation works quite well at low energies, 
but it is not appropriate at higher energies, specially for high Z nuclei|6T].|l5|. Therefore, at higher electron energies 
a different approach is needed to describe the Coulomb distortion effect of the electron. For this purpose, we apply 
the methods of electron scattering where various app roximations have been used to take into account the Coulomb 
distortion effects of the initial and final electron|20(-|2J|. One of them is the Modified Effective Momentum Ap- 
proximation(MEMA) in which the electron momentum and energy are modified by taking into account the Coulomb 
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energy. We have used this approach in the case of charged current quasielastic neutrino scattering and the energy and 
momentum of the electron present in the final state is modified in the Coulomb field of the final nucleus |9j.|l2j-[l4j. 
In the local density approximation, the effective energy of the electron in the Coulomb field of the final nucleus is 
given by: 

E eff = E e + V c (r), 

where 

V c (r) = Z f aAw(- f Pj ^r' 2 dr' + H ^-r'dr') (17) 

Thus, in presence of nuclear medium effects the total cross section a(E v ), with the inclusion of Coulomb distortion 
effects taken into account by Fermi function(MEMA), is written as 

a FF ( MEMA) {Ev) = _ 2G F 2 COS 2 fl c ^ f P ° ^ ^ f 1 



where 



Pe ™« J-i 



-i-VJ^/ F(MEMi) . (18) 



ImU£ F = F{Z, E e )ImU N [E Ve -E e -Q,q\ and 
ImU^ EMA = ImU N [E Uc - E e - Q - V c (r), a] (19) 

III. RESULTS 

In this section we present the results for the total cross section a(E) as a function of energy and the flux averaged 
cross section < a > for various nuclei which have been presently proposed to be studied at SNS facilities using 
neutrinos from stopped muon decays;?]. For the numerical calculation of the cross section, we have classified the 
nuclei in three groups according to their nuclear densities used in this calculation and have been presented in Tables 
I-III. 

In Table-I, we present the nuclear density parameters for 12 C, 14 7V and le O nuclei using a 2-parameter harmonic 
oscillator(H.O.) density given by 

p(r)=p (l + a(-) 2 )exp(-(-) 2 ), (20) 
a a 

In Table-II, we present the nuclear density parameters for 19 F, 23 Na, 27 Al, 28 Si, 31 P, 37 'CI, 40 Ar, 51 V, 52 Cr, 55 Mn, 
56 Fe, 59 C*o, n Ga, 89 Y, 93 Nb, 98 Mo, 115 In, 127 I, 139 La, 181 Ta, 208 Pb and 209 Bi nuclei using a 2-parameter Fermi 
density (2pF) given by 

P(r) = (l + exp((r-a)/a)) (21) 

In Table-III, we present the nuclear density parameters for 32 S, 39 K and 40 Ca nuclei using a three parameter 
Fermi(3pF) density given by 

p(r ) P^ l + W ^) m 

P[ ' ~ {l + exp{{r-a)/a)Y [ ' 

The parameters have been taken from de Vries et al. except for 115 In and 127 1 which have been taken from Ref . |Io| . 
The Q values presented in these tables correspond to the lowest allowed Fermi or Gamow- Teller transitions for the 
above mentioned nuclei except for the case of 40 Ca and 98 Mo for which the Q value corresponding to the ground 
state to ground state transitions have been taken[l6j.|2l|. 
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A. 



Nuclear Medium Effects 



When the reaction v e + n — > e + p takes place in the nucleus, the first consideration is the Q value which inhibits 
the reaction in the nucleus. This inhibition is quite substantial in the low energy region considered here for the nuclei 
like 12 C, 16 <3, 18 Si, 32 S and 40 Ca for which the Q values are rather large (Qsal3-18MeV). In addition to this, the 
effect of Pauli blocking which is taken into account through the imaginary part of the Lindhard function is to further 
reduce the cross section. Finally the renormalisation of weak coupling constants which is generated in our model 
through RPA correlations and is taken into account by calculating the cross section with the modified hadronic tensor 
Jrpa a ^ so re duces the cross sections. In Fig. 4, we have shown the reduction due to these effects separately for some 
representative nuclei like 12 C, 56 Fe, 127 1 and 208 Pb in various mass ranges. 

We see that at low energies the major suppression in the cross section comes due to the consideration of Q-values 
and Pauli blocking in the nuclear medium. The reduction in the cross section a(E) due to these effects decreases with 
the increase of energy. For example at E v — 50MeV, this suppression is ~ 93% for 12 C and « 75 — 77% for other 
nuclei like 56 Fe, 127 1 and 208 Pb respectively (compare the solid lines with the dashed line in Fig. 4). This suppression 
reduces to 40 — 45% in all these nuclei at E v = 200MeV(not shown in Fig. 4). 

In addition to the Pauli blocking, the consideration of RPA correlation in the nuclear medium gives rise to further 
reduction which increases with the mass number and decreases with the increase in energy (compare the dashed 
line with the dotted line in Fig. 4). For example at E v — 50MeV the RPA correlations give a further reduction of 
50% for 12 C, 60% for 56 Fe and around 70% for 127 1 and 208 Pb. As the energy increases it becomes smaller and at 
E v = 200MeF the reduction is 35% for 12 C, 40% for 56 Fe and around 50% for 127 1 and 208 P6(not shown here). 
It should be noted that 40 — 60% reduction due to the medium polarisation effects calculated through the RPA 
correlations in our model is similar to using = 0.7 in some shell model calculations p7) - 130| . 



The effect of Coulomb distortion is calculated using Fermi function F(Z, E e ) as well as with the modified effective 
momentum approximation(MEMA). The results for some representative nuclei like 12 C, 56 Fe, 127 1 and 208 Pb in 
various mass range are shown in Fig. 5. The general effect of the Coulomb distortion of the electron is to increase the 
cross section which depends upon the incident energy and the charge of the final nucleus. 

For a fixed Z, this increase in the cross section decreases with the increase in energy while for a fixed energy the 
inclusion of Coulomb distortion increases with the charge Z. For example for 12 C this is 15% at E v — 50MeV which 
becomes 10% at E v — 200AfeV r . For high Z nuclei the Coulomb effect is very large and results in manifold increase 
in the cross sections. This can be seen by comparing the cross section without Coulomb effect shown by dotted lines 
and the cross sections with Coulomb effects using the Fermi function F{Z 1 E e ) shown by dashed lines. For example 
in the case of 56 Fe nucleus the increase due to Coulomb distortion is 83% at E v = 50MeV which becomes 75% at 
E v = 200A/eV^. However, as discussed in Section-II, the use of Fermi function to calculate the Coulomb distortion 
effects overestimates the cross sections and is not appropriate at higher electron energies. Therefore, we use the 
modified effective momentum approximation(MEMA) and present the results for a(E) in Fig. 5 with solid lines. It 
is seen that for low Z nuclei, like 12 C. the results for the cross sections in the two approximations are quantitative 
similar but the MEMA gives slightly higher cross sections for E v < QOMeV. As the Z increases the cross sections 
with MEMA is higher than Fermi function at lower energies and becomes smaller than the cross sections obtained 
with the Fermi function as the energy becomes large. This energy dependence of the cross section is explicitly shown 
for higher Z nuclei like 56 Fe, 127 1 and 208 Pb in Fig.5. 

We show the results of the energy dependence of the cross sections for all nuclei listed in tables I-III, in Figs. 6-9, 
where dashed lines show the cross sections without RPA correlations and dotted lines show the cross sections with 
RPA correlations using Fermi function for the Coulomb distortion. The solid lines show the cross section with RPA 
correlations where the Coulomb distortion effects are calculated with MEMA. In these figures a comparison of dashed 
lines and dotted lines shows the effect of RPA correlations while a comparison of dotted lines and solid lines shows 
the effect of Coulomb distortion calculated using Fermi function and MEMA. 



B. Effects of Coulomb distortion 



C. Flux averaged cross sections 



We calculate the flux averaged cross section < a > defined as 




(23) 



8 



where (f>rE v ) is given by Eqn.l. The results for < a > are presented in Table-IV, where we show by < a > B j p ^ A the 
flux averaged cross sections with nuclear medium effects without any Coulomb distortion effects. When Coulomb 
distortion effects are taken into account the cross sections without RPA correlations are shown by < a >q and the 
results with RPA correlations are shown by < a > PPA ■ We evaluate < a >q and < a > PPA in a hybrid model where 
at lower energies a{E v ) calculated with the Fermi function and at higher energies a(E u ) calculated with MEMA is 
used to perform the flux averaging in Eqn.23. For low mass nuclei like 12 C, le O, etc. the flux averaged cross section 
has been evaluated with the cross section <j(E) calculated with the Fermi function. Thus, in hybrid model, it is the 
lower value of the cross section which is used for calculating < a >. We see from this table that the effect of the 
Coulomb distortion is to increase the cross section and this increase is quite large for high Z nuclei like 56 Fe, 20S Pb, 
etc. In case of 12 C, it is small but plays an important role in explaining the experimental result (Compare column 1 
and column 3 in table IV) . A comparison of column 2 and column 3 in this table shows the strong reduction due to 
RPA correlations which increases with mass number. 

In Table-V, we compare our results with the results of some other calculations. In this ener gy r e gion of the neutrinos, 
there are many theoretical calculations done for the inclusive neutrino reactions in 12 c@-m|, iHl-lil, while there 
are few calculations for le O 32], 56 Fe 27\, 36] and 20S Pb 27], 30]- 33] and some other nuclei 10]. Some calculations are 
similar to the calculations presented in this pa per Jiot-llll w hile others make use of Shell Model Esl-lpol. r andom phase 
approximation(RPA) with pairing correlations |27j-|28j . |3 1 | - 133 | and elementary particle approach|34j-|36l] . We see that 
for 12 C and 56 Fe our results are in fair agreement with the experimental results and other theoretical calculations. 
For 208 Pb nucleus our results for < a > is comparatively smaller than the results of Refs. [271] . pfflj ] and js^. This 
is mainly due to the different approaches of taking into account the nuclear effects. However, among the different 
calculations of the inclusive cross section < a > in 2t)s Pb, the results do not agree among t hemselves [27j . [30 | - [33 | . 
Therefore, more work is needed for calculating the cross section in 20S Pb at low energies. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the charged current v e reactions on various nuclei which are of present interest. The cross section 
calculations are performed in a local density approximation taking into account the Pauli blocking, Fermi motion 
and the renormalization of weak transition strengths in the nuclear medium. The effect of Coulomb distortion for 
the charged lepton while coming out of the nucleus is taken into account by using the Fermi function as well as the 
modified momentum approximation(MEMA). The cross sections are then averaged over the v e spectra obtained from 
the muons decay at rest where the maximum energy of neutrinos is 52.8MeV. 

We find that 

1. The role of nuclear effects like Q value, Pauli blocking and Fermi motion is to reduce the cross sections. For a 
given Z, this reduction becomes smaller with the increase in energy. For example, at E v — 50MeV, this suppression 
is « 93% for 12 C and w 75 - 77% for other nuclei like 56 Fe, 127 1, etc. This suppression reduces to 40 - 45% in all 
these nuclei at E v = 200MeV. 

2. There is a further reduction of the cross section due to the renormalization of weak transition strengths in the 
nuclear medium. For a given Z, this reduction becomes smaller with the increase in neutrino energy, while for a given 
neutrino energy E v , this reduction increases with Z. For example at E v = 50MeV the RPA correlations give a further 
reduction of 50% for 12 C, 60% for 56 Fe and around 70% for 127 1 and 208 Pb and at E u = 200MeV this reduction is 
35% for 12 C, 40% for m Fe and around 50% for 127 1 and 208 Pb. 

3. The two methods of treating the Coulomb distortion give similar results for low energy neutrinos in the case of 
low mass nuclei. For intermediate and heavy mass nuclei the cross sections with Fermi function are smaller than the 
cross sections with MEMA upto certain energy E Vc after which the the cross sections calculated with Fermi function 
become larger. At the energy E Uc where this cross over takes place changes with nuclei. For example it is around 
40MeV for nuclei like 56 Fe in the intermediate mass range and around 18MeV for nuclei in the heavier mass range 
like 208 Pb. 

4. The total cross sections averaged over the neutrino spectrum obtained from the muons decaying at rest is 
presented for all nuclei considered here. The results for 12 C, 16 0, 56 Fe and 20S Pb nuclei are compared with the 
available experimental results as well as different theoretical calculations. New results have been presented for many 
other nuclei. 

The numerical results presented in this paper can be a very useful benchmark for neutrino nucleus cross section 
measurements being proposed at SNS facilities using various nuclei as nuclear targets. 
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TABLE I: Q-value and the density parameters used in the numerical evaluation of the cross sections, a and a are the parameters 
for HO density used in Eqn.(20). 



Nucleus 


Q (transition) 


a 


a 




(MeV) 




(fm) 


12 C 


17.84 (0+ -> 1+) 


1.69 


1.07 


14 N 


5.60 (1+ -> 0+) 


1.76 


1.23 


16 Q 


19.76 (0+ -> 1+) 


1.80 


1.52 



TABLE II: Q-value and the density parameters used in the numerical evaluation of the cross sections, a and a are the parameters 
used in Eqn.(21) for the 2pF density. 



Nucleus 


Q (transition) 


a 


a 


Nucleus 


Q (transition) 


a 


a 




(MeV) 


(fm) 


(fm) 




(MeV) 


(fm) 


(fm) 


19 F 


3.75(i + -i + ) 


2.59 


0.56 


59 Co 


1.92 (I" _§-) 


4.08 


0.57 


23 N a 


4.56 (f + - | + ) 


2.81 


0.54 


71 Ga 


0.74 (§"->£-) 


4.44 


0.58 


27 Al 


5.32 (§+->§+) 


3.07 


0.52 


89 y 


4.44 (!-->§-) 


4.76 


0.57 


2S Si 


14.80 (0+ -> 0+) 


2.93 


0.57 


93 Nb 


0.91 (| + - 1 + ) 


4.87 


0.57 


31p 


5.91 (| + -i + ) 


3.21 


0.56 


98 Mo 


2.20(0+ -» 6+) 


5.10 


0.57 


37 Cl 


1.32(| + -| + ) 


3.53 


0.52 


115 /n 


0.008 (|+ |+) 


5.36 


0.56 


40 Ar 


4.3 (0+ 1+) 


3.39 


0.61 


127j 


1.42 (f + _ | + ) 


5.40 


0.55 


51 V 


i.26 (r-n 


3.94 


0.50 


139 La 


0.78 (| + - | + ) 


5.71 


0.53 


52 Cr 


5.22 (0+ -f 1+) 


4.01 


0.50 


181 Ta 


0.70(r-|") 


6.38 


0.64 


55 Mn 


0.74 (§--§") 


3.89 


0.57 


208 pb 


5.20 (0+ -f 1+) 


6.62 


0.50 


56 Fe 


6.82 (0+ -f 1+) 


3.97 


0.59 


209 Bi 


3.80(2" -2") 


6.75 


0.47 
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FIG. 4: Quasiclastic cross section a(E v ) ~ £7„ for neutrino reaction in some nuclei. The dashed (dotted) lines represent 
neutrino nucleus cross section without (with) RPA correlations to be compared with the free nucleon cross section (solid lines). 




FIG. 5: Quasiclastic cross section a{E v ) ~ E v for neutrino reaction in some nuclei. The dashed (solid) lines show the cross 
section calculated with RPA correlations and Coulomb effects using Fermi function (MEMA) to be compared with the cross 
section without Coulomb effects (dotted lines). 
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FIG. 6: Quasielastic cross section cr{E v ) ~ E v for neutrino reaction in various nuclei. The dotted (solid) lines show the cross 
section calculated with RPA correlations and Coulomb effect using fermi function (MEMA) to be compared with the cross 
section calculated with Coulomb effect using Fermi function and without RPA correlations (dashed lines-shown here after 
multiplying by 0.6). 
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FIG. 7: Quasielastic cross section a(E v ) ~ E v for neutrino reaction in various nuclei. The dotted (solid) lines show the cross 
section calculated with RPA correlations and Coulomb effect using Fermi function (MEMA) to be compared with the cross 
section calculated with Coulomb effect using Fermi function and without RPA correlations (dashed lines-shown here after 
multiplying by 0.6). 
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FIG. 8: Quasielastic cross section cr{E v ) ~ E v for neutrino reaction in various nuclei. The dotted (solid) lines show the cross 
section calculated with RPA correlations and Coulomb effect using Fermi function (MEMA) to be compared with the cross 
section calculated with Coulomb effect using Fermi function and without RPA correlations (dashed lines-shown here after 
multiplying by 0.6). 
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TABLE III: Q-value and the density parameters used in the numerical evaluation of the cross sections, a, a and w are the 
parameters used in Eqn.(22) for the 3pF density. 



Nucleus 


Q (transition) 


a 


a 


w 




(MeV) 


(fin) 


(fm) 




32 S 


13.20 (0+ -» 1+) 


3.50 


0.63 


-0.25 




7.04 (| + -> | + ) 


3.74 


0.58 


-0.201 


40 C a 


14.80(0+ -» 4") 


3.67 


0.58 


-0.1017 



o.s 




10 20 30 40 10 20 30 40 



FIG. 9: Quasielastic cross section a(E v ) ~ E v for neutrino reaction in various nuclei. The dotted (solid) lines show the cross 
section calculated with RPA correlations and Coulomb effect using Fermi function (MEMA) to be compared with the cross 
section calculated with Coulomb effect using Fermi function and without RPA correlations (dashed lines-shown here after 
multiplying by 0.6). 
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TABLE IV: Total cross section < a >(in 10~ 42 cm 2 ) averaged over the Michel spectrum. < a > I ^f L is the averaged cross 
section calculated with RPA correlations without the Coulomb effect, < a >c PA is the averaged cross section calculated with 
RPA correlations with Coulomb effect and < a >q is the averaged cross section calculated without RPA correlations with the 
Coulomb effect. 



Nucleus 


^ „ \ RPA 
< & >NC 


<o>% 


, „ ^RPA 


Nucleus 


, „ ^ RPA 
< & >NC 


<o>% 


, „ . RPA 


fc 


11.80 


28.10 


13.60 




157.10 


643.18 


254.70 




44.60 


107.91 


52.65 


II Mn 


250.40 


1144.52 


412.40 


fo 


12.00 


30.43 


14.55 


U Fe 


161.40 


685.25 


277.00 


fF 


77.80 


201.07 


96.35 


i?Co 


244.70 


1148.62 


418.70 


llNa 


94.70 


263.50 


120.50 


l\Ga 


335.70 


1762.91 


596.00 


HAl 


82.17 


247.13 


111.62 


89 v 
39 1 


307.10 


1912.11 


633.50 


28 c- 
14 Sl 


33.17 


95.18 


46.25 


%Nb 


370.00 


2332.71 


760.10 


31 p 


89.34 


276.60 


127.10 


llMo 


417.60 


2661.25 


854.00 


32 o 
16 1> 


44.10 


141.38 


64.40 


\fln 


516.60 


3836.37 


1133.00 


nci 
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531.18 


216.00 


127 t 
53 1 


545.40 


4262.30 


1253.00 


18 Ar 


166.50 


560.49 


228.00 


If La 


581.50 


A 7S7 A Q 


i /inn nn 


39 K 
19" 


78.80 


285.14 


123.55 


\fTa 


907.30 


7912.80 


2358.00 


foCa 


38.00 


135.32 


61.25 


IfPb 


902.10 


7857.37 


2643.00 


llv 


198.40 


879.18 


323.20 


IfBi 


824.38 


8257.85 


2497.00 






TABLE V: < a > (10 


42 cm 2 ) for the inclusive reaction for some nuclei. 






Nucleus 




Experimental 
Results 






Theoretical 
Results 




Present 
Calculation. 
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14.40[3JJ, 12.90, 17.60f32] 
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19 Q 










16.90, 17.20[22] 
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